 &MODEL
 NLIST   =           1,
 NSTEP   =        2000,
 EPSI    =  1.0000000E-07,
 WGHT    =   1.000000    ,
 NOBS    =         120,
 NVAXZ   =           4,
 NVAX    =           4,
 NORDER  =           2,
 NVAZ    =           0,
 NLAMZ   =           0,
 NLAMY   =           1,
 NLAMX   =           4,
 NSCALE  =           1,
 NEXPY   =           0,
 NRES    =           0,
 NCOR    =           0,
 NCOV    =           0,
 VARNAME = BC5+A1A2P,
 NVARIANT        =           7
 /
 &AUTOCOR
 RHO     = 2*0.000000000000000E+000  ,
 PI      = 2*0.990000000000000       ,
 NTPI    = 2*1,
 MATFILE = R1.in                                                                           R2.in                                     
                                       
 /
 &YBOXCOX
 VLY     =  0.000000000000000E+000,
 NTLY    =           1,
 NUMLXY  =           0
 /
 &XBOXCOX
 NUMX    =           1,           2,           3,           4, 46*0,
 NUMLX   =           1,           2,           3,           4, 46*0,
 VLX     = 10*0.000000000000000E+000  ,
 NTLX    = 4*1, 6*0
 /
 &YXZDATA
 YXZFILE = Canada1976.in                                                                   ,
 YXZFMT  = (5e15.5)                                                                        
 /

VARIABLE       MINIMUM         MEAN          MAXIMUM   NULL OBS.      SCALE
DEPENDENT
TOTFLOW     0.2577000D+04  0.1476588D+06  0.2736000D+07     0     0.1000000D-05
INDEPENDENT
POP         0.8591300D+05  0.5297039D+06  0.2802800D+07     0     0.1000000D-05
REVE        0.3698300D+04  0.4258776D+04  0.4808200D+04     0     0.1000000D-02
LANGFR      0.5600000D+01  0.7329333D+02  0.9990000D+02     0     0.1000000D+00
UTILITY     0.5655600D-04  0.4915871D-01  0.1172800D+01     0     0.1000000D+01
CONSTANT    0.1000000D+01  0.1000000D+01  0.1000000D+01     0     0.1000000D+01


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) POP         REVE        LANGFR      UTILITY     CONSTANT  
X(l)
===================================================================================================================================

POP          0.1000D+01  0.7451D+02  0.3456D+03  0.4914D+07  0.1889D+06

REVE         0.1342D-01  0.1000D+01  0.4639D+01  0.6595D+05  0.2536D+04

LANGFR       0.2893D-02  0.2156D+00  0.1000D+01  0.1422D+05  0.5466D+03

UTILITY      0.2035D-06  0.1516D-04  0.7034D-04  0.1000D+01  0.3845D-01

CONSTANT     0.5293D-05  0.3944D-03  0.1829D-02  0.2601D+02  0.1000D+01


INITIAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                                                                DERIV.E(Y)      ELAST.E(Y)


POP           X1  0.1389205D+01  0.5457770D-01     25.4537107  0.4002112D-09  0.3872510D+00      1.3892049
                                                                              0.4970914D+00      1.3892049

REVE          X2  0.8322593D+00  0.6932686D+00      1.2004861 -0.6521993D-09  0.2885581D+02      0.8322593
                                                                              0.3626500D-01      0.8322593

LANGFR        X3  0.6643968D-01  0.5550501D-01      1.1970033 -0.8357800D-09  0.1338513D+03      0.0664397
                                                                              0.1916400D-02      0.0664397

UTILITY       X4  0.6335200D+00  0.1496708D-01     42.3275532  0.2744231D-08  0.1902915D+07      0.6335200
                                                                              0.2843971D+02      0.6335200

CONSTANT      X0  0.4955077D+00  0.9808520D+00      0.5051808 -0.4526935D-09  0.7316608D+05      0.4955077
                                                                              0.9221363D-01      0.4955077

VAR(W)            0.1312324D+00  0.1694202D-01      7.7459667  0.0000000D+00

INITIAL VALUES
 
RHO1           0.0000000D+00
RHO2           0.0000000D+00
PI1            0.9900000D+00
PI2            0.9900000D+00
 
LAM(Y)         0.0000000D+00
LAM(X1)        0.0000000D+00
LAM(X2)        0.0000000D+00
LAM(X3)        0.0000000D+00
LAM(X4)        0.0000000D+00


LOG-LIKELIHOOD =  0.3396262D+03

NUMBER OF OBSERVATIONS              =            120

PEARSON-R2(E) (unadjusted)          =         0.8850
PEARSON-R2(E) (adjusted for D.F.)   =         0.8810

PSEUDO-R2(L)  (unadjusted)          =         0.9441
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9422

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.1860993D+00

MEAN of residual W                  = -0.4950736D-12
VARIANCE of residual W              =  0.1312324D+00
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.7100297D+00   Upper Bound =  0.7100297D+00




initialization cards for control routine follow - -  id card - no. parameters , no. steps , step size 
minimization tolerance - starting values of parameters - relative uncertainties of parameters (error estimates) -  
and starting direction  (relative parameter increments for 1st step)
identification line
         9      2000   0.1E+01   0.1E-06
  0.00000E+00  0.00000E+00  0.45951E+01  0.45951E+01  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00

first entry to fcn , parameters of current step follow - - f,x(i)
  -0.3396262428D+03  0.00000D+00  0.00000D+00  0.45951D+01  0.45951D+01  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00
                     0.00000D+00

new absolute minimum , parameters for fcn entry     6  follow - - f,x(i)
  -0.3619399940D+03  0.73463D-02  0.57464D-02  0.45951D+01  0.45951D+01  0.51193D-01  0.97284D-02  0.82294D-05 -0.61471D-03
                    -0.80486D-01

new absolute minimum , parameters for fcn entry    10  follow - - f,x(i)
  -0.3653949242D+03  0.34374D-01  0.28196D-01  0.45951D+01  0.45951D+01  0.96557D-01  0.62246D-01  0.20780D-04 -0.47310D-02
                    -0.61389D-01

new absolute minimum , parameters for fcn entry    14  follow - - f,x(i)
  -0.3720558766D+03  0.11186D+00  0.11219D+00  0.45951D+01  0.45951D+01  0.82560D-01  0.23776D+00  0.70143D-04 -0.19395D-01
                    -0.74724D-01

new absolute minimum , parameters for fcn entry    18  follow - - f,x(i)
  -0.3773067875D+03  0.32090D+00  0.38521D+00  0.45950D+01  0.45953D+01  0.68497D-01  0.20599D+00  0.19904D-03 -0.57616D-01
                    -0.78424D-01

new absolute minimum , parameters for fcn entry    21  follow - - f,x(i)
  -0.3774693713D+03  0.25472D+00  0.44461D+00  0.45927D+01  0.45969D+01  0.75141D-01  0.20957D+00  0.70009D-03 -0.89373D-01
                    -0.75525D-01

new absolute minimum , parameters for fcn entry    25  follow - - f,x(i)
  -0.3784726547D+03  0.16161D+00  0.42348D+00  0.45486D+01  0.46384D+01  0.76836D-01  0.18443D+00  0.10207D-01 -0.95700D+00
                    -0.65069D-01

new absolute minimum , parameters for fcn entry    28  follow - - f,x(i)
  -0.3786732201D+03  0.15461D+00  0.44878D+00  0.45269D+01  0.46691D+01  0.74352D-01  0.18886D+00  0.16135D-01 -0.86966D+00
                    -0.58613D-01

new absolute minimum , parameters for fcn entry    30  follow - - f,x(i)
  -0.3787022504D+03  0.17409D+00  0.45529D+00  0.44961D+01  0.47199D+01  0.74367D-01  0.19394D+00  0.24632D-01 -0.85958D+00
                    -0.58312D-01

new absolute minimum , parameters for fcn entry    32  follow - - f,x(i)
  -0.3787101173D+03  0.16441D+00  0.46767D+00  0.43996D+01  0.48773D+01  0.76377D-01  0.19569D+00  0.52181D-01 -0.88371D+00
                    -0.57266D-01

new absolute minimum , parameters for fcn entry    41  follow - - f,x(i)
  -0.3798868397D+03  0.21716D+00  0.46831D+00 -0.41759D+01  0.17672D+02  0.54248D-01  0.17916D+00  0.23907D+01 -0.11948D+01
                    -0.69346D-01

new absolute minimum , parameters for fcn entry    44  follow - - f,x(i)
  -0.3799489262D+03  0.23004D+00  0.46531D+00 -0.54318D+01  0.19536D+02  0.51408D-01  0.17674D+00  0.27318D+01 -0.12190D+01
                    -0.70860D-01

new absolute minimum , parameters for fcn entry    49  follow - - f,x(i)
  -0.3843607111D+03  0.49192D+00  0.27178D+00 -0.81629D+01  0.24099D+02  0.64213D-01  0.17688D+00  0.35876D+01 -0.83826D-01
                    -0.63232D-01

new absolute minimum , parameters for fcn entry    53  follow - - f,x(i)
  -0.3854326354D+03  0.60462D+00  0.25214D+00 -0.96773D+01  0.26606D+02  0.70573D-01  0.18262D+00  0.40580D+01 -0.37554D+00
                    -0.62124D-01

new absolute minimum , parameters for fcn entry    57  follow - - f,x(i)
  -0.3855765767D+03  0.65525D+00  0.25494D+00 -0.10112D+02  0.27327D+02  0.67564D-01  0.18151D+00  0.41949D+01 -0.24140D+00
                    -0.64705D-01

new absolute minimum , parameters for fcn entry    60  follow - - f,x(i)
  -0.3855839278D+03  0.65354D+00  0.25467D+00 -0.10131D+02  0.27349D+02  0.66305D-01  0.18279D+00  0.42026D+01 -0.25040D+00
                    -0.66338D-01

new absolute minimum , parameters for fcn entry    62  follow - - f,x(i)
  -0.3855847642D+03  0.65337D+00  0.25597D+00 -0.10163D+02  0.27398D+02  0.66409D-01  0.18170D+00  0.42166D+01 -0.24809D+00
                    -0.66768D-01

new absolute minimum , parameters for fcn entry    64  follow - - f,x(i)
  -0.3855848846D+03  0.65413D+00  0.25474D+00 -0.10207D+02  0.27466D+02  0.66443D-01  0.18172D+00  0.42341D+01 -0.24576D+00
                    -0.66843D-01

new absolute minimum , parameters for fcn entry    67  follow - - f,x(i)
  -0.3855857859D+03  0.65392D+00  0.25619D+00 -0.10732D+02  0.28269D+02  0.66857D-01  0.18293D+00  0.44540D+01 -0.24500D+00
                    -0.66707D-01

new absolute minimum , parameters for fcn entry    71  follow - - f,x(i)
  -0.3855905070D+03  0.65438D+00  0.25463D+00 -0.13841D+02  0.33018D+02  0.66238D-01  0.18081D+00  0.57645D+01 -0.25987D+00
                    -0.66600D-01

new absolute minimum , parameters for fcn entry    73  follow - - f,x(i)
  -0.3855917951D+03  0.65293D+00  0.25579D+00 -0.15190D+02  0.35078D+02  0.66411D-01  0.18166D+00  0.63361D+01 -0.24859D+00
                    -0.66675D-01

new absolute minimum , parameters for fcn entry    75  follow - - f,x(i)
  -0.3855918747D+03  0.65333D+00  0.25564D+00 -0.15680D+02  0.35827D+02  0.66443D-01  0.18158D+00  0.65422D+01 -0.25116D+00
                    -0.66657D-01

new absolute minimum , parameters for fcn entry    77  follow - - f,x(i)
  -0.3855918752D+03  0.65332D+00  0.25563D+00 -0.15703D+02  0.35862D+02  0.66431D-01  0.18156D+00  0.65521D+01 -0.25105D+00
                    -0.66661D-01

new absolute minimum , parameters for fcn entry    78  follow - - f,x(i)
  -0.3855918752D+03  0.65331D+00  0.25563D+00 -0.15703D+02  0.35862D+02  0.66431D-01  0.18156D+00  0.65521D+01 -0.25105D+00
                    -0.66660D-01

new absolute minimum , parameters for fcn entry    79  follow - - f,x(i)
  -0.3855918752D+03  0.65331D+00  0.25563D+00 -0.15703D+02  0.35862D+02  0.66431D-01  0.18156D+00  0.65521D+01 -0.25105D+00
                    -0.66660D-01


the function (fcn) has not improved by    0.10000D-06 after 3 steps which include a minimum - - -


therefore the run has been terminated
  -0.3855918752D+03  0.65331D+00  0.25563D+00 -0.15703D+02  0.35862D+02  0.66431D-01  0.18156D+00  0.65521D+01 -0.25105D+00
                    -0.66660D-01


all numbers of last entry refer to final pt. of run


the lowest value of the function (fcn) was  -0.3855918752D+03 for entry   79 - - x(i) follows
  0.65331E+00  0.25563E+00 -0.15703E+02  0.35862E+02  0.66431E-01  0.18156E+00  0.65521E+01 -0.25105E+00 -0.66660E-01



run statistics follow - -

 nfcn=   83 nst=   25


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) POP         REVE        LANGFR      UTILITY     CONSTANT  
X(l)
===================================================================================================================================

POP          0.1000D+01  0.2332D+02  0.8525D+03  0.3713D+07  0.7512D+06

REVE         0.4289D-01  0.1000D+01  0.3656D+02  0.1593D+06  0.3222D+05

LANGFR       0.1173D-02  0.2735D-01  0.1000D+01  0.4355D+04  0.8812D+03

UTILITY      0.2693D-06  0.6279D-05  0.2296D-03  0.1000D+01  0.2023D+00

CONSTANT     0.1331D-05  0.3104D-04  0.1135D-02  0.4943D+01  0.1000D+01


FINAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T      STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                   (CONDIT.)        (CONDIT.)      (CONDIT.)                   DERIV.E(Y)      ELAST.E(Y)
                                                    (PAR = 0)      (PAR = 1)

POP           X1  0.2881043D+00  0.2756455D+00      1.0451985                -0.2966108D-09  0.3986930D+00      1.4302514
                                 0.9599019D-02     30.0139311                                0.3655287D+00      1.4984851

REVE          X2  0.9819764D-24  0.3822723D-21      0.0025688                -0.7976113D+12  0.9295850D+01      0.2681109
                                 0.1688779D-23      0.5814712                                0.9058621D+01      0.3075248

LANGFR        X3  0.1093417D+01  0.1974937D+01      0.5536465                -0.1521741D-10  0.3398986D+03      0.1687153
                                 0.1840948D+00      5.9394255                                0.4735142D+03      0.1951207

UTILITY       X4  0.8889819D+00  0.4031481D+00      2.2051001                 0.4895647D-10  0.1480375D+07      0.4928477
                                 0.3893679D-01     22.8314119                                0.3589627D+08      0.6737598

CONSTANT      X0  0.4472426D+01  0.5779499D+01      0.7738433                -0.5807204D-11  0.2995049D+06      2.0283579
                                 0.1010201D+01      4.4272648                                0.2751954D+06      2.2131519

VAR(W)            0.1880258D+00  0.1238067D+00      1.5187055                 0.0000000D+00
                                 0.2853406D-01      6.5895223
 
RHO1      (E)     0.5738953D+00  0.1105361D+00      5.1919267                 0.3951232D-04
                                 0.1003518D+00      5.7188330
RHO2      (E)     0.2502072D+00  0.1436671D+00      1.7415759                 0.3349563D-04
                                 0.1362786D+00      1.8359982
PI1       (E)     0.1514013D-06  0.1491000D+00      0.0000010     -6.7069083 -0.3784969D+02
                                 0.1390918D+00      0.0000011     -7.1894951
PI2       (E)     0.1000000D+01  0.1172220D+01      0.8530824      0.0000000  0.5113949D+01
                                 0.1130522D+01      0.8845470      0.0000000
 
LAM(Y)    (E)     0.6643091D-01  0.3214172D-01      2.0668124    -29.0453966 -0.1026197D-03
LAM(X1)   (E)     0.1815610D+00  0.7293478D-01      2.4893612    -11.2215179  0.2473408D-04
LAM(X2)   (E)     0.6552148D+01  0.4639164D+02      0.1412356      0.1196799 -0.1356795D-04
LAM(X3)   (E)    -0.2510545D+00  0.5907373D+00     -0.4249850     -2.1177847  0.4260142D-05
LAM(X4)   (E)    -0.6666040D-01  0.2420076D-01     -2.7544753    -44.0754912  0.4016481D-04


LOG-LIKELIHOOD = -0.1272269D+04

NUMBER OF OBSERVATIONS              =            120

PEARSON-R2(E) (unadjusted)          =         0.9953
PEARSON-R2(E) (adjusted for D.F.)   =         0.9948

PSEUDO-R2(L)  (unadjusted)          =         0.9829
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9808

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.1466759D+06

MEAN of residual W                  = -0.4408935D-01
VARIANCE of residual W              =  0.1880258D+00
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.8939835D+00   Upper Bound =  0.8058048D+00


CPU TIME:      3.39 sec
